Tunneling into Nonequilibrium Luttinger Liquid with Impurity 
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We evaluate tunneling rates into/from a voltage biased quantum wire containing weak backscat- 
tering defect. Interacting electrons in such a wire form a true nonequilibrium state of the Luttinger 
liquid (LL). This state is created due to inelastic electron backscattering leading to the emission 
of nonequilibrium plasmons with typical frequency hui < U. The tunneling rates are split into two 
edges. The tunneling exponent at the Fermi edge is positive and equals that of the equilibrium LL, 
while the exponent at the side edge Ef — U is negative if Coulomb interaction is not too strong. 

PACS numbers: 73.23.-b, 71.10.Pm, 73.21.Hb, 73.40.Gk 



By virtue of advances in modern nanotechnology elec- 
tron tunneling spectroscopy became a powerful tech- 
nique that enables to reveal electron correlations on a 
mesoscale. Suppression or enhancement of the tunneling 
conductance at low bias is a signature of electron inter- 
action in the system and is commonly called a "zero-bias 
anomaly" (ZBA) [J,i]. Measurements of ZBA in disor- 
dered metals in high-mobility two-dimensional elec- 
tron gases Q , in the edges of quantum Hall systems [j| , 
and recent measurements of magneto-tunneling in arrays 
of quantum wires Q are milestones of this field. 

Of particular interest on this way is the study of elec- 
tron tunneling into quantum nanowires 0,0,1. Cen- 
tral to much of the fascinating physics of these one- 
dimensional (ID) electron systems is that Coulomb in- 
teraction has a dramatic effect leading to the emergence 
of the Luttinger liquid (LL) 9]. This strongly corre- 
lated state of matter is commonly described in terms of 
bosonic elementary excitations. Measurements of ZBA's 
in nanowires confirm predictions based on the LL model. 

Notably, the behavior of a strongly correlated quantum 
system can change drastically when it is driven out of 
equilibrium. Remarkable examples include the Kondo 
phenome na flC L 11 , 12, IH and the Fermi-edge singularity 
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FIG. 1: Tunneling experiment with a voltage biased quantum 
wire (x > 0). Right moving electrons have a larger chemical 
potential relative to left moving electrons, hr — = U > 0. 



elastic electron backscattering at the defect induces the 
emission of real nonequilibrium plasmons with typical fre- 
quencies fko < U. In the non-dissipative LL they transfer 
the shot noise of backscattered current with the Poisso- 
nian statistics to the distant point of tunneling x (Fig. [I]), 
thereby considerably influencing the ZBA. Let us empha- 
size an important difference between the present setup 
and that of Ref. 0. While in the model of [12] the 
nonequilibrium state is "injected" into the LL, here it is 
created by a scatterer located inside the LL. 

Our results can be summarized as follows. We consider 
a spinless LL with e-e interaction described by the con- 
ventional parameter K. The tunneling rates into/from 
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problem [U,|15(. Recent experiments initiate the study thfi right electron states ( R _ stat es) in the LL are split 



of the nonequilibrium tunneling spectro scop y of carbon 
nanotubes [lg] and quantum Hall edges [17| . 

In this paper we consider the tunneling into a volt- 
age biased one-channel ballistic wire, containing a weak 
backscattering defect (Fig. 1). Previous studies of this 
model focused on the nonlinear conductance and shot 
noise [2l| . However, the tunnel spectroscopy 

of this problem, which requires the analysis of the single- 
particle Green's function, has never been addressed. We 
show that interacting electrons in such a wire form a 
generic nonequilibrium LL state, characterized by non- 
Gaussian plasmon correlations, and develop a real-time 
instanton approach to evaluate the tunneling rates. In- 



into two edges, r|(e) = r±(e) + Ar ± (e), as shown 
in Fig. [2] The first term here accounts for the equi- 
librium contribution to ZBA around the Fermi energy, 
r^(e) x6(±e)\e/E F \ 21 , with exponent 7= {l-K) 2 /AK. 
The nonequilibrium corrections are singular at e = — U: 
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in case of tunneling from the R-state into the tip and 
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in case of tunneling from the tip into the R-state. Here 
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FIG. 2: Energy dependence of electron tunneling rates from 
the nonequilibrium LL (left pane) and into the LL (right 
pane), shown for rfj = 0.2 and different strengths of repulsive 
e-e interaction: (1) K = 0.4; (2) K = 0.5; and (3) K = 0.75. 



5 = (1 - K)/2 and r% = r 2 (E F /U) 4S (T(2K) < 1 is 
the renormalized reflection coefficient 23] . The tunneling 
rates into/from the left electron branch remain almost 
equal to the equilibrium ones, T^(e) ~ Y^(e + U). The 
above power-law singularities are smeared on a scale of 
the nonequilibrium dcphasing rate 



= (2/7r)J7r^sin 2 ttS. 
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It is worth stressing the oscillatory dependence of 
on the interaction parameter K. It differs from that ob- 
tained in the model of Ref . [22[ , which reflects a different 
type of the nonequilibrium LL state. 

The result ([2|) corresponds to inelastic tunneling with 
absorption of real plasmons. An electron tunneling into 
the LL with the energy e < can accommodate itself 
above the Fermi energy of right moving states by picking 
up the energy ftw > \e\ from the nonequilibrium plasmon 
bath. Since the energy of out-of-equilibrium plasmons 
is limited by the applied voltage, one has a threshold: 
e > —U. The correction fll} describes the inverse pro- 
cesses — inelastic tunneling from the LL accompanied 
by the stimulated emission of nonequilibrium plasmons 
with typical energy fko ~ U. In the absence of interaction 
the splitting of tunneling rates T^(e) can be understood 
as the result of the double-step distribution of R-states 
due to the scattering off the impurity [241 ]. 

We model the wire as the spinless LL [9] with lin- 
ear dispersion and short-range forward e-e interaction, 
characterized by the amplitude Vq, which fixes K — 
(1 + Vq/ttvf)^ 1 / 2 . We decouple the electron field into 
right- and left-moving fields ip±. The wire contains a 
weak impurity at x — with a bare reflection amplitude 
r <C 1. 

Our theoretical analysis is based on the functional 
bosonization (FB) [2^]. We consider an electron motion 
in a Hubbard-Stratonovich field ip which mediates e-e in- 
teraction. The special feature of ID geometry is that by 
a local gauge transformation, ip v i— > e ii/j^, the coupling 



between ip± and ip can be removed everywhere except 
the points of scattering provided that d^fi^ — —tp, where 
d n = d t + t]Vf9 x . The phases $± define charge and cur- 
rent responses to the potential *p in each chiral branch, 
Pr] = d x -d v /2ir and j v = —d t "d v /2TT. Following this gauge 
transformation one constructs the bosonized Keldysh ac- 
tion S = Sb + <Simp in terms of the variables 1? and <fi I2II 



S b 



/ dt dxY 









2 \ ttvf J 



(4) 



Simp = --|r | 2 Tr 



- „*4 - — i* 

9+e g~e 



x=0 



Here v = (2-kvf)^ 1 is the ID density of states, $ = 
d- — and g± are "quasiclassical" Green's function in 
the right/left leads. They are fixed by boundary condi- 
tions g n = (1 - 2f v )(h + ih) - h, f±(h - h) being the 
electron distribution functions. In the limit of zero tem- 
perature one has f±(t) = ie~ v± */27r(t + iO), where the 
chemical potentials satisfy /x+ — fj,- = U. The diagonal 
matrix $ = diag($~, $ + ) and the Pauli matrices f act 
in the Keldysh space, the upper indices ± referring to 
the two branches of the Keldysh contour C. The trace 
operation (Tr) is performed in the Keldysh x time space. 
The quadratic action Sb describes the charge and current 
fluctuations in the clean wire, while the impurity action 
Simp accounts for plasmon emission and absorption due 
to electron backscattering in the lowest order in ro| 2 . 

We start by considering Gaussian fluctuations of 9± 
and ip, described by the action Sb- Within the FB elec- 
tron phases have no free dynamics — in contrast to the 
conventional bosonization — but rather respond to the 
internal electric field. This response is found by opti- 
mizing Sb for a given tp, which gives the gauge rela- 
tion drj'drj = —(p. One has to solve it by taking the 
proper structure of the Keldysh theory into account, 
$T)[<p] — —D 0ri T 3 ip. We have introduced doublets, e.g. 
= (<p~ ,(p + ) T , and the bare particle-hole propagator 
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E F 1 is a short time 



cut-off. Then the quadratic action Sb, expressed solely in 
terms of <p, assumes the RPA form, Sb = ^<p T V~ 1 (p : with 
a nonlocal effective interaction V = Vq — (Vq /ir)K 2 d x V. 
Here V = IXLa/^-VA) an d is the propagator of the 
plasmon modes moving with velocity u = vp/K, 
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where c+ = 1 
tin 



■7, c 



= 7 and C ±_ Al = (l-^ 2 )/(4^). 
To find the tunneling rates we represent the electron 
Green's function at the point of tunneling x > as a 
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path integral over the field tp, 

G^(x,r) = f ®tpe iS ^ iS ^G}{x,T,Q; [tp]). (7) 
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Here G v (x, t, t'; [tp]) denotes the Green's function for 
a given configuration of tp. It satisfies the Dyson 
equation with the spatially local self-energy T* v [<p] = 
— iS(x)(\ro\ 2 vp/2) e 171 ® e~"'* due to impurity scatter- 
ing, the phase [ip]) being the linear functional of tp 
introduced above. The action 5j describes the creation 
of a hole at time t = and an electron at the instant 
t = r, while Sf corresponds to the inverse process, 



Sf{T,X,[<p])=tf(x,T)-0*(x,O) 



dtdx^J^.(8) 



The source here, e.g. J> = 8(x — x) (S(t), —S(t — r)) , 
acts on both branches of C as shown in Fig. [3] 

To find the Green's function we proceed with a semi- 
classical approximation 27]. One looks for a saddle- 
point trajectory tp* which optimizes the total action 
<Stot = S + Sj and further estimates the tunneling rate by 
evaluating <Stot [<P*]- For <C 1 we can find such a trajec- 
tory approximately imposing that it minimizes only the 
quadratic part of the action, Sq = Sb + Sj, which gives a 
simple linear equation in tp* . Taking into account correc- 
tions to tp* of order of |r| 2 , which follow from the exact 
non-linear equations of motion, would lead to a contribu- 
tion ~ |r| 4 to the tunneling action, which is beyond the 
accuracy of our method. 

Using the above approximation we find ^[y*] = 
3> m J. Here the phase-phase correlation function satisfies 
the relations v d t = Do^^x-V^ and #^(0, 0) = 0, 
that enables easy evaluation of •&* using the Eqs. (0) 
and |6]). For instance, in case of tunneling from the 
tip into the right branch (77 = +) the relative phases 
<J>* = d*_ — 7?+ explicitly read 

t + x/u =p ia 



»$?(*) =ln 



t — t + x/u =F ia 
S 
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(9) 



t — t — x/u -f- i(x\ /t — r + x/vp^fia 
t — x/u — ia 



t + x/vf =f ia 
into the RPA action we ob- 



Substituting tp* 

tain iSQ-(r) = — 27 ln(±ir/a). The impurity action 
evaluated on the instanton ([9]) consists of four terms, 
S* mp ( T ) = Ea/3 5 im P ( r )' the indices a,(3 = ± arising 



from the Keldysh structure of 
contribution is given by 



At U > the main 
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Here we have modified the l/£ 2 -behavior of the bare equal 
point polarization operator by taking into account the 
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FIG. 3: Visualization of Sj [ Eq.© ] on the Keldysh contour. 



time-dependent LL renormalization of the reflection am- 
plitude r(t) ~ ro(t/a) 2S , which results from quantum 
fluctuations around the saddle-point trajectory tp*. 
We further concentrate on the limit of long x, so that 
<C x/vp — x/u, where = max{r0,l/Ae} is the 
typical accommodation time and Ae = min{|e|, [U + e|} 
is the energy relative to the nearest edge. In this case 
the instanton ([9]) consists of well separated plasmon and 
particle- hole kinks, moving with velocities u and vp, and 
gives two independent contributions, S* and S* h , to the 
impurity action. Their long-time behavior at t > 1/U is 
defined by singularities of the integrand (fTO]) : 

iS;(r) ~ rid (iUr) 25 e WT -\r\/2 Tli> , (11) 
iS* ch {r) ~ r 2 u {~ l Ur) iS e-^T{2K), 

where the rate t^ 1 is given by Eq. ([3]) and the numerical 
factor Ci = T{2K)/T 2 (\ + 5). This asymptotics is iden- 
tical in both cases of tunneling into and from the LL. 
The linear growth of S* (r) stems from the time domain 
1 1 1^2 1 Tj while the oscillations cx T 2S e lllT are governed 
by the far distant times |ti 2 + x/u\ < r. The action 
S* reveals the Poissonian statistics of the shot-noise of 
backscattered current carried by plasmons (Fig. [4j). Sim- 
ilar to that the action S* h describes the shot-noise due 
to inelastically excited electron-hole pairs at the distant 
times \t1_2 + x/vp\ < r. 

Within the saddle-point approximation we obtain from 
the representation 

G%{x,t) ~ -iwjIfTje^W^^W^W. (12) 

Tunneling rates are related to G^{x,t) by the Fourier 
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FIG. 4: Real plasmons (wavy line) created in the course of 
inelastic electron backscattering at t ~ —x/u, are absorbed 
by the injected electron at t ~ 0; typical duration of a scat- 
tering/tunneling event is At ~ r <C x/u. 
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FIG. 5: Differential tunnel conductance between the BCS- 
type tip and the nonequilibrium LL. The parameters are: U = 
A, rl = 0.2, (I) K = 0.4; (2) K = 0.5; and (3) K = 0.75. 



transformation. Since <C 1, one can expand the os- 
cillatory part of the action in Eq. (|12|) , retaining only 
the first term. We have checked that the particle- 
hole contribution 5* h exactly cancels the lst-order im- 
the Green's function, AG + = 
v , where the average is performed 



purity correction to 
(G , + £ + (cp)G , + e iS J) 
with the RP A- action 5b. Keeping then only the plasmon 
contribution to the tunneling action, we finally find 
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{U/E F f-'T{-2 1 ) Im{(Tz) 



2~, 



(13) 



-C 1 (±I) 2 ''*(-2 7 , 1 + 25- 27, -I - z)} , 



where z = (e + i/2r^)/U is the complex energy and \t 
is the confluent hypergeometric function [2{| . The latter 
is singular at z — > —I, yielding the power laws stated in 
Eqs. (U) and ©. We plot the rates (H2J) in Fig. [2] versus 
energy e for different strengths of e-e interaction, indi- 
cating the edge exponents. Remarkably, in the vicinity 
of the edge e = — U the in-rate is enhanced, while the 
out-rate is suppressed, provided the nonequilibrium 
exponent A = 2(7 — 5) is negative, which is the case of 
not too strong interaction realized at K > i. 

We illustrate our theory by considering tunneling in 
the LL from a superconducting tip with the singular BCS 
density of states ^t(e) oc |e|(e 2 — A 2 ) -1 / 2 . Current mea- 
surements using this setup enable to reveal the nonequi- 
librium structure in the tunneling rates [l6j |. For the 
tunneling current we have 



I = Gn 



de 



fs-vTt -(l-h-v)T- itie-V), (14) 



where Tf = T^(e) + L^(e) and Gt is the bare tunnel 
conductance, fn Fig. [5] we show the differential con- 
ductance dl/dV in units of the normal state conduc- 
tance Gt(A) ~ Gt(A/^) 27 at the scale A. Due to 
the double-edge structure of the tunneling rates (fl~3|) the 
peaks of the BCS density of states are split by the bias 
voltage U and show power-law behavior with exponents 



Ai = 27 - 1/2 and A 2 = 2(7 - 6) - 1/2. Singularities 
of u t (e) visibly enhance the nonequilibrium structures in 
the rates Tf, making the conductance profile strongly 
asymmetric, even in the limit of small rfj. 

To summarize we have developed a real-time instanton 
approach to the problem of tunneling into the nonequi- 
librium state of the interacting quantum wire containing 
weak backscattering defect. Tunneling rates are split into 
two edges, the power-law exponent A at the nonequilib- 
rium edge e = — U being negative, provided the repulsive 
e-e interaction is not too strong (K > i). This nonequi- 
librium effect is associated with inelastic electron tunnel- 
ing accompanied by absorption/emission of real plasmons 
with a typical frequency Hlu ~ U. The approach devel- 
oped in this work will be useful for analysis of tunneling 
and interference in a broad class of nonequilibrium LL 
structures with impurities and/or tunneling couplings. 
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